Pure spin current generation in a Rashba-Dresselhaus quantum channel 
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We demonstrate a spin pump to generate pure spin current of tunable intensity and polarization in 
the absence of charge current. The pumping functionality is achieved by means of an ac gate voltage 
that modulates the Rashba constant dynamically in a local region of a quantum channel with both 
static Rashba and Dresselhaus spin-orbit interactions. Spin-resolved Floquet scattering matrix is 
calculated to analyze the whole scattering process. Pumped spin current can be divided into spin- 
preserved transmission and spin-flip reflection parts. These two terms have opposite polarization of 
spin current and are competing with each other. Our proposed spin-based device can be utilized for 
non-magnetic control of spin flow by tuning the ac gate voltage and the driving frequency. 

PACS numbers: 73.23.-b, 73.21.Hb, 72.25.Dc, 72.30.+q 
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I. INTRODUCTION 

Manipulation of electron spins can be achieved via 
applying external active control, which is the essen- 
tial requirement of spintronics devices X Especially, spin- 
resolved current generation is one of the key interests in 
spintronics research for its potential application in quan- 
tum information science? 2 "! 3 . Various approaches were pro- 
posed to overcome the fundamental challenge in the is- 
sues of spin current manipulation, detection, and injec- 
tion efficiency. Methods based on controlling magnetic 
fieldi£i& and material ferromagnetism 7 - are investigated. 
However, for practical applications, more efficient meth- 
ods that do not involve strong magnetic field or inter- 
faces between ferromagnets and semiconductors are still 
needed. Spin pumping can be a viable solution to the 
spin current generation*^ 

Pumping of charge current is a fully quantum mechan- 
ical phenomenon in a mesoscopic system that can gener- 
ate current without applied bias between two leads. The- 
oretically and experimentally, charge current pump has 
been realized and implemented in a quantum channel or 
a cavity in the way of periodic modulation ! 10 ' 11 ! 12 : 13 ! 14 
In the adiabatic regime, Brouwer proposed a clear pic- 
ture that the pumped current depends on the enclosed 
area in parametric space which is formed by a set of pe- 
riodically varied parameters. Such formalism was read- 
ily extended to non-adiabatic regime, which is valid in 
the whole spectrum of frequency^ If the spin degree 
of freedom is incorporated, spin-dependent transmission 
coefficients can be differentiated either directly by ex- 
ternal magnetic fiel d 16 : 1 or by spin-orbit interaction^^ 
Spin pumping is generalized from the quantum pump- 
ing and exempted from the spin injection problem which 
occurs in the integrated semiconductor-fcrromagnet ar- 
chitecture. 

In order to achieve spin pumping, a Rashba-type nar- 
row channel (which ignores the presence of the Dressel- 
haus term) driven by local time-dependent potential was 
proposed ] 19 : 20 When electrons propagate through the po- 
tential region, quasi-bound state feature was shown to 



enhance the spin-resolved transmission difference so that 
sizable pure spin current can be generated. However, 
since the Dresselhaus spin-orbit interaction is an intrin- 
sic effect in semiconductor materials with bulk inversion 
asymmetry— it is essential to take into account this effort 
when considering such a spin pumping device. It should 
be noted that the presence of the Dresselhaus term will 
lead to the spin-flip mechanism which can modify the 
spin-pumping characteristics in a qualitative way. We 
shall elucidate the possibility to manipulate not only the 
intensity but also the polarization of the spin current. 

In this paper, the spin-resolved Floquet scattering ma- 
trix formalism is applied to our system i 22 : 23 Based on 
the Floquet theorem, this formalism provides an ex- 
act and nonperturbative solution to the time-periodic 
Schrodinger equation in the mesoscopic system. Be- 
cause the time-dependent spin-orbit interaction couples 
two spin polarizations and all sidebands together, ana- 
lytic expression for the sideband dispersion is not feasible. 
Thus, we determine the sideband dispersion relation nu- 
merically by solving the Schrodinger equation in a nearly 
complete basis. Besides, the spatial inhomogeneity can 
also be handled by matching boundary conditions piece 
by piece spatially. The Floquet scattering matrix gives a 
coherent solution that goes beyond the adiabatic regime. 



II. MODEL AND FORMALISM 

The system under consideration is a two-dimensional 
electron gas (2DEG) that is present at the interface of 
a heterostructure due to modulation doping and has in- 
trinsic static Rashba and Dresselhaus spin-orbit interac- 
tions. The system configuration is shown in Fig. [T] A 
quasi-one-dimensional (Q1D) narrow channel is formed 
from the 2DEG via a lateral confining potential (along 
the y direction). The barrier separating the Q1D channel 
from the 2DEG should be strong enough so the tunneling 
time between them is much longer than the carrier trans- 
port time in the Q1D channel. A finger gate is placed in 
the middle of the channel (the grey region in Fig. 1) 
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FIG. 1: Schematic illustration of the quasi-lD spin-orbit 
quantum channel embedded in 2DEG. In this narrow channel, 
the electron gas has static Rashba and Dresselhaus spin-orbit 
interactions which are characterized by Qo and /3q respectively. 
The central grey region, with width I, is biased by ac gate volt- 
age so that Rashba strength locally modulated as ai cos(u)t). 
The origin of :r-axis is set at the left edge of the grey region. 



that modulates the local Rashba interaction strength si- 
nusoidally via an ac-bias. Hence, the system can be de- 
scribed by the effective Hamiltonian 
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+ W so (r, t) + V c {y), 



(1) 



where m* denotes the electron effective mass and V c (y) 
indicates the confinement potential in transverse (y) di- 
rection. 7ig* atIC and 7i so (r,i) characterize, respectively, 
the static and dynamic parts of spin-orbit interaction. If 
we consider a narrow quantum channel where the sub- 
band energy spacing is large enough to decouple p y from 
spin-orbit interaction, the intersubband mixing is thus 
neglectedi 19 ' 20 The longitudinal part of the dimension- 
less Hamiltonian is then given by 

Hi = k 2 x - a Q a v k x + l3 cr x k x , 

H x (t) = ~a 1 a v cos(u;t){k x ,e(l/2-\x-l/2\)},(2) 

where a l (i = {x,y,z}) denotes Pauli matrices and k x 
indicates the momentum operator —id x . Anticommuta- 
tor {• ■ ■ } is used to maintain the hermitianity of H. x (t). 
The static Rashba strength ao is proportional the electric 
field perpendicular to the interface where 2DEG lies. Ad- 
ditionally, (3q is the phenomenological Dresselhaus cou- 
pling parameter. In the finger gate region, the Rashba 
parameter oscillates sinusoidally with amplitude ot\. For 
simplicity, we restrict the subsequent discussions to the 
lowest subband and ignore the subband index. The con- 
tributions from other subbands can be added if a more 
realistic consideration is needed. 

To proceed, it is convenient to rotate the spin quantiza- 
tion axis such that 7i°. is diagonalized. The transformed 
Hamiltonian is 



k l - lo° z k x , 
1 



(3) 



cos(ujt){k x , 9(1/2 - \x - l/2\)}, (4) 



where = (a z sin </> — a y cos 0) , 70 = y/otQ + > an d 
<j> — arctan(ao//3o)- W x illustrates not only our choice of 
spin-up and spin-down states but also that the location of 
subband bottom is at — 7q/4. Based on Floquet theorem, 
the wave functions in lead L (x < 0) and lead R (x > I) 
are given by 
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where Xa denotes the spinor basis and fj, represents the 
incident energy. The sideband index m runs essentially 
for all integers. From the dispe rsion relation in Eq. 
©,fc^ lo . and k^ a are |[t? ct 7o + V^o + + muJ )} and 
K^o-To — \/7o + + muj)] respectively, where r\ a is de- 
fined as a* a a j . a TO ,<r (a! m a ) is the amplitude of the right- 
ward (leftward) wave in the rath sideband with spin a in 
lead L. Similarly, b m _ a (b' m <7 ) is for lead R. Technically, 
these amplitudes are determined by boundary condition 
and the direction of incident wave. 

In the time-dependent region M (0 < x < I), the gen- 
eral solution would be 



i(e-\-mi>j)t 



Xo 



(6) 



where e is the Floquet quasi-energy. ^ m ,a-{x) is solved 
from Schrodinger's equation, 



a.\k x 
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a' 

(*m+l,a' + *m-l,<7')] = HHV' ( 7 ) 

These coupled equations can be expressed in matrix form, 
fc^ + ^HWf^H' '*, (8) 

where 



^■(m,ff)(m',o-') — 2 a {v,<T')^ m ' m ' + 1 fim,m'-l) 
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(9) 

(10) 

(11) 



Because this is a transport problem, we have to solve the 
eigenvalue q for fixed e. This quadratic eigenproblem can 
be solved by introducing another of auxiliary equation 
= qty. Then Eq. © becomes 



1 

H (°) -hH I*' 



(12) 



If we truncate the sideband index m at —M/2 and M/2, 
where M is an even integer, the eigenvalues q 3 and eigen- 
vectors "0m a are numerically determined from the above 
secular equation. 
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Because Hamiltonian in Eq. |J5J| preserves time-reversal 
symmetry, any q> is associated with —(q J )*, i.e. e(q^) = 
e{ — (q J )*). In addition, for Hamiltonian is also invariant 
under inversion followed by spin flip, q J has its another 
counterpart — q 1 . Thus, we can definitely sort the (AM + 
4) complex eigenvalues into two groups. 

For the case of evanescent modes, those right-decaying 
waves are characterized by positive Im(q , - y ); left-decaying 
waves have negative Im(g J ). On the other hand, for the 
case of propagating modes that have real g J , we sort qi 
with positive (negative) group velocity to be rightward 
(leftward) propagating waves. The group velocity is de- 
termined by (de/dqi)^ Therefore, the wave function in 
region M is given by 

*W(aO = ' ; ' r + g' j ^e iqi,L % (13) 



where superscripts R and L are added to indicate the 
propagating or decaying direction. 

Wave functions are matched in the time domain by e 
= n and continuous across the boundaries. Their deriva- 
tives satisfy the following boundary conditions: 

a x *(x,t)| a=a f-a ar *(a ! ,t)| 1B= o- - ^coB(wt)<7**(0,i), 

d x t>(x,t)\ x=b - -d x V(x,t)\ x=b+ = ^icos(wt)(7*¥(M). 

(14) 

The above boundary conditions can be written down in 
matrix form, 



a + a' 

(K R a + K L a') - (S R Q R g + S L Q L g') 



S R e 4Q 'g + S L e iQ 



g = 



S R g + S L g', 
i(H( 1 )- 70 S)(a + a'), 



(K R b + K L b') - (S R Q L e l Q R; g + S l QV qL V) = i(H« - 7o S)(e 4Kllz b + e^b'), 



(15) 
(16) 
(17) 
(18) 



where those column vectors a, g, and b, are assigned 
values from amplitudes a m ^, 9j, and 6 mjCT respectively. 
The above (2M + 2) x (2M + 2) matrices' S R(L) , Q R(L \ 
X, and K R ( L ) have matrix elements 



S R ( L ) _ jlJML) 



Q R(L) 

^3,3 



3 ML) 



^(m,cr)(m',ff') — 3m,m'ficr,cr'(~ <^ct,T + ^<r,l)> 

(19) 

After some algebra, we have the following matrix equa- 
tion from Eqs. (T5j) to (fT8|) : 



Mu M12 
M 2 i M22 



(20) 



M 



Mil Mil) denotes ( 4M + 4 ) x (4Af + 4) matrix 
connecting the input coefficients with output coefficients 
including all propagating and evanescent Floquet side- 
bands. 

In order to construct the Floquet scattering matrix, 
we need to introduce the concept of probability flux am- 
plitude into M. We can straightforward define a new 
matrix as 

. (21) 



( V R 



Mu M 12 \ [V R 
M21 M22 { V L 



where V R ^ L \, , „ 



> III .111 U (T,(T 



2k 



R(L) 



In both leads, V R ( L ) takes the form of diagonal matrix 
with the square root of group velocity absolute value from 
each sideband and spin type. It is worth mention that 
M! is not unitary yet due to the presence of evanescent 
modes. In the final stage, we obtain a unitary Floquet 
scattering matrix by setting the evanescent modes of the 
total scattering matrix M! to be zero: 



S 



11 V 
T 1V 



(22) 



The unitarity of Floquet scattering matrix reflects the 
current conservation law , 25 ' 26 and is used as the criteria 
to check numerical convergence. 

The reflection and transmission coefficients are readily 
obtained by summing over matrix elements of S. When 
electrons that are incident from L lead with initial spin 
Oi are partially reflected and transmitted to final spin o/, 
the spin-resolved reflection and transmission coefficients 
are written as 



R afai( £ ) = ^2 l^(m,<r / )(0,ff < )| 
m 

T a f L *i( £ ) = ^ \^m,<rf)(0,<T t )\' 



(23) 
(24) 



On the contrary, if the electron is incident from lead R, 
this gives rise to such reflection and transmission coeffi- 
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cients 



m 

T^(e) = ^|r' (m ^ )(0 , CTi) | 2 . (26) 



Under zero longitudinal bias, the spin-resolved current 
pumped out through lead R is generally defined as 



if — 



def(e)[Tf L + T, 



■-n 



nRR 

TT 



R 



RR 

U 



I J dsf{s)[T™ + T^ + Rf« + Rfx -I], 



(27) 

where /(e) is Fermi-Dirac distribution. The spin-resolved 
current can be derived based on the framework of 
Biittiker's formula 27 by regarding two spin types as dif- 
ferent terminal channels. The generalization of Biittiker's 
formula for Floquet scattering matrix has been strictly 
proven . 28 ' 29 The spin current and charge current at lead 



if and if = If + If. 



R are defined as if 

Because system Hamiltonian in Eq. ^ has inversion fol 
lowed by spin flip symmetry, we can transform the trans- 
mission and reflection coefficients as T^ R . = T^ f _ a , 
and Rff a . — i?^^ / _ cr . . Such transformation firstly guar- 
anteed that there is zero charge current in this system. 
Secondly, when certain amount of spin current is pumped 
out at lead R, there should be equal amount of spin cur- 
rent with opposite polarization pumped out at lead L. 
Furthermore, if this symmetry is combined with current 
conservation condition, spin current formula can be sim- 
plified to a more convenient form in calculation: 

lf = \\ dsf(s)[(Tf L - Tjf ) + (Rff Rf*)}. (28) 

The first two terms represent contributions from trans- 
mitted electrons, and the last two terms are attributed 
to reflected electrons whose spin is changed. Hence, we 
separate if into spin-preserved transmission and spin- 
flip reflection parts because their effects are different and 
discussed in the following context. Thus, 



jR jR, trans 



rii,ren 



(29) 



It should be noted that if there is no Dresselhaus term, 
the lf' reR term is identically zero, and the if is then 
reduced to the same form in Ref. 1 19t 




FIG. 2: Spin-resolved transmission and reflection coefficients 
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R" l , and Tfi as functions of the incident energy. 
The values of reflection coefficients are multiplied by 50 to 
clarify the shape of the curves. Qo = 0.12, /3o = oio, I — 30, 
w = 0.002, and ai =(a) 0.04, (b) 0.06, (c) 0.08. The spin 
current, if 1 , which depend on Fermi energy for different ai 
are plotted in (d). 




FIG. 3: (a) Pumped spin current 7^* the Fermi energy (related 
to the bottom of the first subband in Q1D channel). /3o is 0, 
0.03, 0.06, and 0.12. Parameters ai = 0.08, a = 0.12, I = 30, 
and u = 0.002. (b)/f, if ' trans , and if ' rcfl are illustrated for 
/3q = 0.12 case in (a). 



III. RESULTS AND DISCUSSION 

Utilizing the above derived formula in previous sec- 
tion, it is easy to calculate the spin current pumped from 
the spin-orbit quantum channel via numerical means. 
The reasonable material parameters are chosen from 
the narrow-gap heterostructure based on InGaAs-InAlAs 



based system. According the experimental data, we as- 
sume that the 2DEG has an electron density n e = 1 x 10 12 
cm~ 2 , effective mass m* — 0.04to , and a = 0.12 
(ha = 2.8 x KT 11 eV m)^ The ratio between Rashba 
and Dresselhaus terms can vary in certain range due to 
experimental difficulties. 31 Thus, we examine the cases 
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for flof&o varying between and 1. In our calculations, 
the length and energy units are chosen to be I* = 4.0 nm 
and E* = 59 meV (the Fermi energy of the 2DEG). We 
assume that the ac-biased gate has a width of / = 30/* 
and its driving frequency is chosen as Tllu = 0.002_E* 
(w/2tt — 28 GHz). The bottom of the lowest energy 
level (first subband) in the Q1D channel is assumed to 
be slightly below the Fermi level, E* of the 2DEG so 
that the Fermi energy relative to the bottom of the first 
subband in the Q1D channel (denoted Ep) is compara- 
ble to hto. All numerical results are obtained for zero 
temperature. 

The dependence of transmission and reflection coef- 
ficients on the incident electron energy (e) for various 
values of a% are illustrated in Figs. [Ha)-(c) when the 
static Rashba and Dresselhaus constants are the same, 
i.e. /3o = olq. In order to clarify the important features 
shown in these figures, we redefine the energy zero at the 
bottom of the first subband with the presence of Rashba 
and Dresselhaus terms, i.e. H'° — > H'° +7q/4. The coef- 
ficients T^ L , Tf L L , Rff-, and Rff, which are needed for 

calculating if, are plotted in Figs. 2(a)-(c). For trans- 
mission coefficients, we find sharp features at integer val- 
ues of e/oj, indicative of the resonant inelastic scattering. 
As ai increases, the dip around e/oj = 1 moves toward 
lower energy, and the dip width is broadened. The reason 
for the shift of dip location is that a stronger oscillating 
potential would lower the real part of the quasi-bound 
state energy and shorten the lifetime of electrons trapped 
in such a stated When a\ is increased to 0.08 as shown 
in Fig. 2(c), a higher order resonance seen as a shallow 
dip around e/oj = 2 becomes more apparent because of 
the absorption and emission of two quanta (with energy 
2Tiuo). The most significant effect of the Dresselhaus in- 
teraction is the emergence of the spin-flip process, which 
leads to appreciable spin-flip reflection coefficients, R R ^ 
and Rff. In Figs. 2(a)-2(c), Rff and Rff have a saw- 
like behavior with peaks appearing at integer values of 
e/u>, where electrons are bounced back due to the pres- 
ence of quasi-bound states. Although their values are still 
minute compared with T^ L and T RL , they can lead to 
significant change in the final spin current when we take 
differences of the spin-up and spin-down contributions. 

Figure 2(d) illustrates the spin current as a function of 
the Fermi energy, e p (which reflects the carrier density in 
the Q1D channel) for various values of ai. The curves in 
this Figs 2(a)-(c) can be approximately divided into two 
parts: the low energy region (s/uj < 2) and high energy 
region (e/oj > 2). In the low energy region, the reflection 
coefficients are too small compared to (Tll L — Tf^ ) , and 
If is dominated by the contribution due to transmission 
process (denoted /^. trans ). i n the high energy region, 
{Rf* - Rff ) becomes stronger than (T^ L - Tff) and 
the contribution to if due to reflection process (denoted 
J^ ,refl ) becomes dominant. As ai increases from 0.04 
to 0.06, more spin current is pumped out the first peak 
at Ef/ld = 1 and into the second peak at Ef/lu = 2. 



When a± is tuned even higher to 0.08, high order reso- 
nances become more relevant. Thus we have a further 
enhanced peak around Ef/oj = 2 and a reduced peak 
around Ef/oj = 1- However, because (Rff — Rff') is 
always negative, lf' lcfi results in negative contribution 
to if and it pulls the spin current curves downward. For 
Ep/ui > 2, lf' lefi becomes dominant so that negative spin 
current is generated. As otx increases, if (for Ep/u> > 2) 
becomes more negative due to higher probability of the 
spin-flip process. 

In Fig. 3(a), we focus on the effect of Dresselhaus inter- 
action on the pumped spin current for a fixed ct\. In the 
case of zero 0o, only one kind of spin polarization can be 
pumped. 19 As /3q increases, If curves tend to shift down- 
ward due to increased spin-flip scattering process. In the 
low density case (ef/oj < 2), experimentally reasonable 
Po may hardly change the sign of if. In the higher den- 
sity case (sp /u> > 2), the sign of if is more vulnerable to 
the strength of the Dresselhaus term. When (3 is 0.03, 
0.06, and 0.12, the threshold values of £ at which the sign 
of if starts to change are at ef/u = 4.89, 3.09, and 2.32, 
respectively. 

A simple physical picture is presented here to give a 
conclusive explanation. The conditions in Fig. 3(b) are 
taken as an example. Based on the dispersion relation of 
H'® in Eq. (O, when the electron is incident from lead 
L, \k R ^\ is always larger than il f° r the same energy. 
Thus, it is easier for spin-up electron to tunnel through 
this oscillating barrier due its larger flux, i.e. this dis- 
persion of static Hamiltonian tends to favor T RL rather 
than Tf L . On the other hand, because scattering poten- 
tial can be approximately regarded as proportional 
to momentum, spin-up electrons could be more suscepti- 
ble to the scattering process so that T RL is favored here. 
In low energy region, these two mechanisms are compet- 
ing so that (T& L — T RL ) may be positive or negative and 
jr, trans nag odious peaks. 

In high energy region, because the second mechanism 
is less relevant, only monotonically increasing /^< trans is 
present. For /^' ren ; the situation is just on the opposite 
side. Because \Jsq A is greater than {k^J, there would 
be less chance for incident spin-down electrons to be re- 
flected. Hence, lf' refl always contributes to negative spin 
current and is monotonically decreasing. When incident 
energy is low, lf' lcfi only compensates part of /^ trans . 
When energy increases, lf' rcii becomes dominant and 
there is a threshold e/oj beyond which if starts to change 
sign. 



IV. CONCLUSION 

We have proposed a promising approach to generate 
spin current non-magnetically in the absence of charge 
current. A quasi- ID channel with static Rashba and 
Dresselhaus spin orbit interaction is studied. Spin pump- 
ing is achieved by an ac gate voltage to locally modulate 
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the Rashba constant. Pumped spin current can be at- 
tributed to both the spin-preserved transmission and the 
spin-flip reflection processes. These two terms contribute 
to opposite polarization of the spin current. 

It is found that in the low density case {sf/lo < 2), the 
spin-preserved transmission is dominant and featured by 
resonant inelastic scattering. In the high density case 
{ef/^J > 2), there is a threshold beyond which spin cur- 
rent begins to switch polarization. Furthermore, it is 
found that the static Dresselhaus coefficient (3 as well 
as the dynamic Rashba coefficient ct\ can enhance the 
spin- flip process and modify the threshold value of Ef/uj, 
at which the spin polarization switches. In conclusion, 



we have demonstrated a feasible way to control dynam- 
ically the intensity and polarization of the spin current 
via changing the strength of the ac-biased gate voltage 
and tuning the driving frequency. 
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